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1. Introduction 

Let srf = {Ai, . . . ,A r } be a set of real m x m matrices. As usual, for n > 1, let 
us denote by &/ n the set of all n-products of matrices from ; £/° = I. For each 
n > 1, define the quantity 

/0n GO = max p{A in ■ ■ ■ A i2 A h ), 

where maximum is taken over all possible products of n matrices from the set 
and p(-) denotes the spectral radius of a matrix, that is the maximal magnitude of 
its eigenvalues. The limit 



lim sup {p n {s/)) 



l/n 



is called the generalized spectral radius of the matrix set stf 0, [l(| • 
Similarly, given a norm || • || in M m , for each n > 1, define the quantity 

Pn{#t) = max \\A in ■■■Ai 2 A il \\, 

where \\A\\, for a matrix A, is the matrix norm generated by the vector norm 
in R m , that is \\A\\ = supi| x i| =1 Then the limit 



p(£?) = lim sup (p r 



l/n 
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does not depend on the choice of the norm || • || and is called the joint spectral 



radius of the matrix set stf [291 ] . 



For matrix sets srf consisting of a finite amount of matrices, as is our case, the 
quantities p(^) and p(&/) coincide with each other [5|] and their common value is 
denoted as 

p[sf) = p{${) = p{s/), 

while the quantities p n (^) and p n (£^) form lower and upper bounds, respectively, 
for the joint /generalized spectral radius: 

Pn(*f) < As*) = f>W) < Pn{*f), Vn>0. 

This last formula may serve as a basis for a posteriori estimating the accuracy 
of computation of p(&&)- The first algorithms of a kind in the context of control 
theory problems have been suggested in 0], for linear inclusions in [lj, and for 



problems of wavelet theory in [8|, |9j, Later the computational efficiency of these 
algorithms was essentially improved in 13|, . Unfortunately, the common feature 
of all such algorithms is that they do not provide any bounds for the number of 
computational steps required to get desired accuracy of approximation of p{&&)- 

Some works suggest different formulas to compute p{&^)- So, in [7| it is shown 
that 

p(£&) = lim sup max |tr(Aj n ■ • ■ Ai 2 Ai 1 )\ 1 ' n , 

where, as usual, tr(-) denotes the trace of a matrix. 
Given a norm || • || in M. m , denote 

= max ||>1||. 

Then the spectral radius of the matrix set s$ can be defined by the equality 

9 {fif) = ipfKH, (l) 

where infimum is taken over all norms in R m 12, 2^]. For irreducible matrix setsjl] 
the infimum in ([1]) is attained, and for such matrix sets there are norms || • || in 
M m , called extremal norms, for which 

11*11 < P(*0- (2) 

In analysis of the joint spectral radius ideas suggested by N.E. Barabanov 
play an important role. These ideas have got further development in a variety of 
publications among which we would like to distinguish [3ll |. 

Theorem 1.1 (N.E. Barabanov). Let the matrix set s$ = {Ai, . . . , A r } be irre- 
ducible. Then the quantity p is the joint (generalized) spectral radius of the set srf 
iff there is a norm \\ ■ \\ in M m such that 

p\\x\\ = max \\Aix\\. (3) 



X A matrix set srf is called irreducible, if the matrices from si have no common invariant sub-spaces except 
{0} and R m . In 18-20] such a matrix set was called quasi-controllable. 



Throughout the paper a norm satisfying (J3]) will be called a Barabanov norm 
corresponding to the matrix set s# '. Note that Barabanov norms are not unique. 

Similarly, [26l . Thm 3.3], [28l ] the value of p equals to p(s^) if and only if for some 
central-symmetric convex bodjQ S the following equality holds 

pS = conv(\jA i s) , (4) 



where conv(-) denotes the convex hull of a set. As is noted in [26], the relation 
(j3]) was proved by A.N. Dranishnikov and S.V. Konyagin, so it is natural to call 
the central-symmetric set S the Dranishnikov-Konyagin-Protasov set. The set S 
can be treated as the unit ball of some norm || • || in M. m (recently this norm is 
usually called the Protasov norm). As Barabanov norms as Protasov norms are 



the extremal norms, that is they satisfy the inequality ([2]). In [24 . \2B . \32\ it is 
shown that Barabanov and Protasov norms are dual to each other. 

Remark that formulas ([2]), (J3|) and ([4]) define the joint or generalized spectral 
radius for a matrix set in an apparently computationally nonconstructive manner. 
In spite of that, namely such formulas underlie quite a number of theoretical con- 
structions (see, e.g., [4 lit! 17, 23|, 31, and algorithms [27J for computation of 



Different approaches for constructing Barabanov norms to analyze properties of 



the joint (generalized) spectral radius are discussed, e.g., in [14j,ll5[ and [30, Section 
6.6]. 

In [2~ll ] the so-called max-relaxation algorithm was proposed for computation of 
the joint spectral radius of matrix sets. In the paper an alternative iteration proce- 
dure, a linear relaxation procedure, is introduced that allows to build numerically 
Barabanov norms for the irreducible matrix sets and simultaneously to compute 
the joint spectral radius of these sets. 

The paper organized as follows. In Introduction we give basic definitions and 
present the motivation of the work. In Section [2] the iteration procedures is intro- 
duced. This procedure is called the linear relaxation procedure since in it the next 
approximation to the Barabanov norm is constructed as the linear combination of 
the current approximation and some auxiliary norm. Section [3] is devoted to the 
proof of convergence of the iteration procedure. In Section U] we briefly describe 
the so-called max-relaxation iteration scheme for computation of the joint spectral 
radius. At last, in concluding Section [5] we present results of numerical tests and 
discuss some shortcomings of the proposed approach. 



2. Linear relaxation iteration scheme 

Let g/ = {Ai, . . . ,A r } be an irreducible set of real m x m matrices, || • ||o be a 
norm in M m , and e 7^ be an arbitrary element from lR m satisfying ||e||o = 1. 
Let A - and A + be fixed but otherwise arbitrary numbers satisfying the condition 

< A~ < A + < 1. 

These numbers will play the role of boundaries for parameters of the linear relax- 
ation scheme below. Define recursively the sequence of the norms || • || n , n = 1, 2, . . ., 
according to the following rules: 



The set is called body if it contains at least one interior point. 



LRi: if the noTTfi || * \\ n has been already defined compute the Quantities 

I IXLclXj \\fi . IXLciXj jtL^£C Ij-j ma m / r~\ 

p£ = max 77— r , p n =mm — n , j n = max |^e|| R ; (5) 

x '7^0 \\ x \\n x ¥=0 \\ x \\n » 

LR2: choose an arbitrary number A n € [A~, A + ] and define the norm \\ ■ || n +i- 

||x||„+i = A„||x|| n + (1 - A n )7^ 1 max ||^ix|| n . (6) 

The iteration procedure ([5]), © will be referred to as the linear relaxation pro- 
cedure {the LR-procedure) since in it the next approximation ||x|| n +i to the Bara- 
banov norm is constructed as the linear combination of the current approximation 
1 1 a; I | n and some auxiliary norm. 

As we will see in Section 13.11 p~ < p < p+ for any n = 0, 1, . . ., and so the 
quantities {p~ } form lower bounds for the joint spectral radius p of the matrix set 
8$ ', while the quantities {p„} form upper bounds for p. 

Remark that the norm ([6]) is correctly defined for any choice of j n because due 
to irreducibility of the matrix set srf = {^4i, . . . ,A r } for any 1 / the vectors 
Aix, . . . , A r x do not vanish simultaneously, and then p~ > as well as 7„ > 

Pn \\ e \\n > 0. 

Before we start proving that the LR-procedure converges to some Barabanov 
norm and that the quantities p^ converge to the joint spectral radius p of the 
matrix set srf make two remarks. 

Remark 1. The norms || • \\ n satisfy the normalization conditions ||6||n — 1? ^ — 
1,2, ... , which can be derived by induction from ([6]). Then by © 

max, LAje L 
In = n— n 



and therefore 

In e [p~,pt], n = 0, 1,... . (7) 

Remark 2. Instead of the iteration procedure ([5]), ([6]) one can consider the fol- 
lowing, formally more general, procedure in which the quantities 7„ are chosen 
arbitrarily if only they satisfy the inclusions (J7|), and the obtained norms are nor- 
malized forcibly: 

LR' X : provided that the norm || • 1 1 ^ has been already found compute the Quantities 

1 maxj ||Ajx|| n _ . max, ||j4jx|| n . . 

pi = max — n , p„ = mm — n (8) 

x^O \\x\\ n x^O \\x\\ n 

LRg: choose arbitrary numbers A n G [A~,A + ], 7„ £ [PniPn] an d build first the 
auxiliary norm \\ ■ ||° +1 : 

IMIn+l = ^rall^lln + (1 — A n )7 n 1 max ||j4jx|| n , 

i 

and then define the norm \\ ■ \\ n +i in such a way that the normalization condition 
IMU+i = 1 be satisfied: 

IMU+l = IMIn+l/IHIn+l- ( 9 ) 



I 



In fact, if to write down formulas for recalculation of the norms ||x|| n +i via 
and to represent them in the form similar to ©: 

IMU+i = A^||x|| n + (1 - A^)(7^) _1 max ||A;x|| n , 

% 

then one can find that the corresponding quantities will be uniformly separated 
from zero and unity while the quantity j' n will be equal to the quantity "f n defined 
by ©• The corresponding calculations are not complicated but cumbersome and 
because are omitted. 

So, consideration of the iteration procedures of the form (|HJ), © gives nothing 
new, and such procedures are not studied in what follows. 



3. Proof of the main result 

Clearly, to prove that the iteration procedure ([5]), © converges to some Barabanov 
norm || • ||* (and that the quantities p^ converge to the joint spectral radius p of 
the matrix set &/) it suffices to prove Assertions Al, A2 and A3: 

Al: the sequences {p^} an d {Pn} are convegent; 

A2: the limits of the sequences {p^} and coincide: 

p = lim pX = lim p~; 

n—toc n— >oo 

A3: the norms \\ ■ \\ n converge pointwise to a limit \\ ■ \\* . 

Properties of the iteration procedure (jSJ), ((6]) needed to prove Assertions Al, A2 
and A3 are established below. 



3.1 Relations between and p 

Lemma 3.1. Let a,/3 be numbers such that in some norm \\ ■ \\ the inequalities 

a\\x\\ < max ll-Aixll < /Sllxll, 

hold. Then a < p < (3, where p is the joint spectral radius of the matrix set s/ . 

Proof . Let || • ||* be some Barabanov norm for the matrix set Since all norms 
in R m are equivalent, there are constants o~ > and a + < oo such that 

o-~\\x\\* < \\x\\ < a + \\x\\*. (10) 

Consider for each k = 1,2,... the functions 

A k (x) = max \\A ik . . . A^A^xW- 

l<ii,i 2 ,---,ifc<'' 

Then, as is easy to see, 

a k \\x\\ < A k {x) < /3 k \\x\\. (11) 



rr 



Similarly, consider for each k = 1,2,... the functions 



A fcO) == . max \\A ik . . . A i2 A ilX \\*. 

I<ii,i2,...,ik<r 

For these functions, by definition of Barabanov norms the following identity hold 

Al(x) = p k \\x\\*, (12) 



which is stronger than (jlip . 

Now, note that (jlOp and the definition of the functions A) t (x) and A£(x) imply 



a-At(x) < A k (x) < a + AUx). 



Then, by ([II]), (HI, 



a V < p k < — f3 k , Vfc, 



from which the required estimates a < p < f3 follow. □ 

So, Lemma [37X1 and the definition ([5]) of imply that the quantities {Pn} form 
the family of lower bounds for the joint spectral radius p of the matrix set g/, while 
the quantities {p^} form the family of upper bounds for p. This allows to estimate 
a posteriori errors of computation of the joint spectral radius with the help of the 
iteration procedure ©-(H]). 



3.2 Convergence of the sequence of norms {|| • \\ n } 

Given a pair of norms || • ||' and || • ||" in R m define the quantities 

e-(|M|',|M|'') = minW;, e+(|| . ||',|| . \\>) = ™ M. ( 13 ) 

Since all norms in M m are equivalent to each other, the quantities e~(|| • ||', || • ||") 
and e + (|| • ||', || • ||") are correctly defined and 

o<e-(ii-ir,ii-ir')<e + (ii-ir,n-ir')<oo. 

Therefore the quantity 

eocfli • w, ii • n = ^14445 > i, (14) 



which is called the eccentricity of the norm || • ||' with respect to the norm || • ||" 
(see, e.g., (12]), is also correctly defined. 
Let us start proving convergence of the sequence of the norms || • || n . 

Lemma 3.2. Let \\ ■ \\* be a Barabanov norm for the matrix set si . Then the 
sequence of the numbers ecc(|| • || n , || • ||*) is nonincreasing. 

Proof. Denote by p the joint spectral radius of the matrix set s/. Then by defini- 
tions of the function e + (-) and of the Barabanov norm || • ||* from the relations ([5]), 



([6]) we obtain: 



|x|| n+ i = A„||x|| n + (1 - \n)l n 1 max \\AiX\\ n < 



from which 



|n+l, 



|*) (A n ||x|r + (1 - A n )7" 1 max H^zf) = 

= • ||n, II • ID (AnIMP + (1 " Xn)ln VlMI*) , 

*)<e + (||-|| n ,||-|r)(A n + (l-A n ) 7 -V). (15) 



Similarly, by definitions of the function e (•) and of the Barabanov norm || • ||* 
from the relations ([5]), ([6]) we obtain: 

||x||„+l = A n ||x|| n + (1 - Xnhn 1 max \\AiX\\ n > 

i 

> e '(\\ ■ lln, II • II*) (An||x||* + (1 - A n ) 7 " 1 max||Aix|r) = 

= e ~(ll • lln, || • ID (AnlMI* + (1 - Khn 1 P\M*) , 

from which 

e"([| • ||n + i, || • ID > e-(|| • || n , || • |D (A n + (1 - A n ) 7 - 1 p) ■ (16) 
By dividing termwise the inequality (|15p on (|16p we get 

/ii II || n*\ e+ (ll ' l|n+l, || ' || ) . e+ (ll ' lln, || ' || ) /ii n n |i*n 

ecc( " " 11 ' 11 } = Firram - Frem = ecc( » • ii • ii } - 

Hence, the sequence {ecc(|| • || n , || • ||*)} is nonincreasing. □ 

Denote by A r i oc (R m ) the topological space of norms in R m with the topology of 
uniform convergence on bounded subsets of W n . 

Corollary 3.3. The sequence of TIOTTYIS {|| • I |n} ^ COTYipCLCt ifl A^ oc (IR m ). 

Proof. For each n and any x ^ by the definition (I13j) of the functions e + (-) and 
e~(-) the following relations hold 

e-(||-||n,||-|r)<|{^<e + (||-||n,||-|r), 

and then 

e-(\\-LA\-\n<¥k<e + (\\-\u\-n 

||e|| 

from which 

ecc(|| - |L, II - II*) llelh " e "™ " " X " n " eCC( " ' " ' "^Lf" 6 ""- 



Since here the norms || • ||^ by Remark [Tl satisfy the normalization condition ||6| 
1, and by Lemma l3T2l ecc(|| • || n , || • ||*) < ecc(|| • ||o, || • ||*), we obtain 

1 ll-^ll * ll^ll * 

II lux OU ^ W X Wn < ecc(|| • ||o, || • H* ) 17-777 • 



ecc(J| • Ho, || • \D || „ 

Therefore the norms || • || n , n > 1, are equicontinuous and uniformly bounded on 
each bounded subset of M m . Moreover, their values are also uniformly separated 
from zero on each bounded subset of M. m separated from zero. From here by the 
Arzela-Ascoli theorem the statement of the corollary follows. □ 

Corollary 3.4. If at least one of subsequences of norms from {|| • || n } converges 
in iVi oc (IR m ) to some Barabanov norm then the whole sequence {||*||n} also converges 
in iVi oc (lR m ) to the same Barabanov norm. 

Proof . Let {|| • ||„ fe } be a subsequence of {|| • || n } which converges in iVi oc (lR m ) to 
some Barabanov norm || • ||*. Then by definition of the eccentricity of one norm 
with respect to another 

ecc (ll ' \\n k -, || • ID — > 1 as k — > oo. 

Here by Lemma 13.21 the eccentricities ecc(|| • || n , || • ||*) are nonincreasing in n, and 
then the following stronger relation holds 

ecc( || • \\ n , || • ||*) — > 1 as n — > oo. (17) 

Note now that by the definition (|13p . (|14p of the eccentricity of one norm with 
respect to another 

' < < eccfll • ILJ • n, 



ecc(|| • || n , || • ||*) 

from which by ()17f) it follows that the sequence of norms {|| • || n } converges in space 
Ni oc (W m ) to the norm || • ||*. □ 

Lemma 3.5. Assertion A3 is a corollary of Assertions Al and A2. 

Proof . By Corollary 13.31 the sequence of norms {|| • || n } has a subsequence {|| • ||„ fe } 
that converges in space A r i oc (lR m ) to some norm || • ||*. Then, passing to the limit 
in (J3J) as n = — > oo, we get by Assertions Al and A2: 

maxj \\Aix\\* 

P= Tj— 77; , 

||x||* 

which means that || • ||* is a Barabanov norm for the matrix set £/. This and 
Corollary 13.41 then imply that the sequence {|| • || n } converges in space A r i oc (IR m ) to 
the Barabanov norm || • ||*. Assertion A3 is proved. □ 

In view of Lemma [3. 51 to prove that the iteration procedure ©, © is convergent 
it suffices to verify only that Assertions Al and A2 hold. 



3.3 Convergence of the sequences {p^} 

In the same way as in Section HI from Lemma 13.11 and the definition ([5]) of it 
follows that quantities {p^} form the family of lower bounds for the joint spectral 



o 



radius p of the matrix set srf , while the quantities {Pn} form the family of upper 
bounds for p. This allows to estimate a posteriori errors of computation of the joint 
spectral radius with the help of the iteration procedure dSJ) , dSJ) - 

To prove that the sequences {p^} are convergent, let us obtain first some auxiliary 
estimates for maxj ||^4jx|| n+ i. By definition, 



max ||Ajx|| n+ i = max \ A n ||^4ix|| n + (1 - X n )^ n 1 max || AjAix\\ n } . 
i * [ 3 J 



(18) 



Here for each i the summand (1 — \n)l n 1 max^ ||j4jj4jx|| n in the right-hand part is 
estimated, by the definition ([5]) of the quantities p^, as follows: 

p~(l - A n )7~ 1 ||Ax||„ < (1 - Khn 1 max\\AjAix\\ n < p+(l - X n h~ l \\Aix\\ n . 

3 

Therefore 



max{A n ||Aix|| n + p n (l- A n )7 T , 1 ||Aiar|| Jl } < 
i 

< max {\ n \\Aix\\ n + (1 - A n )7~ 1 max ||^Ajx|| n | < 

< max{A n ||Aix|| n + p+(l - A n )7 r ^ 1 ||^ i x|| n } . (19) 

i 

Here by the definitions (J3J), © of the quantities p~ and of the norm ||x|| n +i we 
have 

max{A n Pix|| n + p~(l - A n )7^ 1 ||A i x|| n } = 

i 

= (A n + p~(l- Xnhn 1 ) max||Ajx|| n = 
= A n max ||Aix|| n + p~(l - Xnhn 1 max ll^lln > 

i i 

> p~A n ||x|| n +p~(l ~ A n )7^ 1 max|| J 4 i x|| n = p~||x|| n+ i. (20) 

i 

Similarly, by the definitions ([5|), © of the quantities /?+ and of the norm ||x|| n +i 
we have 

max {\ n \\Aix\\ n + pt(l ~ Kh^W^Wn} = 

i 

= (A n + p+(l - Xnhn 1 ) max \\AiX\\ n = 

= X n max ||Aix|| n + p+(l - A,^" 1 max ||^4ix||„ < 

i i 

< p+A n ||x|| n + p+(l - Xnhn 1 max \\Aix\\ n = p+||x|| n+ i. (21) 



The estimates (fT8j) - (|2Tp imply 

p~||x|| n+ i < max Pjx|| n+ i < /9+||x|| n+ i, 



n 



from which 



and then 



- , maxj \\AiX \\ n +l , + w , n 

Pn < inj <Pn. V x ¥= 0, 

Flln+1 



P n < P n+1 < P n+l < Pn- 



So, the following lemma is proved. 

LEMMA 3.6. The sequence {p^} is bounded from above by each member of the 
sequence {p^} and is nondecreasing. The sequence {p^} is bounded from below by 
each member of the sequence {p^} and is nonincreasing. 

In view of Lemma 13.61 there are the limits 

p~ = lim p~, p + = lim p+ 

n— Yoo n— voo 

which means that Assertion Al holds. Hence, to prove that the iteration procedure 
©, © is convergent it remains only to justify Assertion A2: p~ = p + ■ 

To prove that p~ = p + below it will be supposed the contrary, which will lead 
us to a contradiction. 



3-4 Transition to a new sequence of norms 

To simplify further reasoning we will switch over to a new sequence of norms for 
which the quantities p^ will be independent of n. 

As was established in Corollary 13.31 the sequence of the norms || • || n is compact in 
space N\ oc (W n ). Consequently, there is a subsequence of indices {n^} such that the 
norms || • \\ nk converge to some norm || • ||* satisfying the normalization condition 
||e||o = 1 while the the quantities X Uk and -y nk converge to some numbers po and 
770 respectively. Then, passing to the limit in ([5]), by Lemma 13.61 we obtain: 

I IIIQjX i A. i x\\r\ niclX^ n ITlclX^ -A-i 6 n 

/r = max r—r- , p = mm n — — , rjo = n~~ rrs • 

x^o \\ x \\o \\ x \\o \\ e \\o 

Now by induction the following statement can be easily proved. 

Lemma 3.7. For each n = 0, 1, 2, . . . the sequence of the norms \\ ■ \\n k +n converges 
to some norm \\ ■ ||* satisfying ||e||* = 1, and the sequences of the quantities X nk+n 
and 7n fe +n converge to some numbers p n 6 [A~, A + ] and rj n respectively. Moreover, 
for each n = 0, 1, 2, . . . we have the equalities 

maxj LAjX * ■ . max, LAjX * _ maxj LA^e * / n n\ 
m l x flu =p ' m l n ij — ij« =p ' iTih =r?n ' ( ' 

and the recurrent relations 

IMIn+l = MnlMIn + (1 ~ Pn)"n~ X max (23) 



Note that the norms (|23[) and Q are correctly defined since, by irreducibility of 
the matrix set stf = {A\, . . . , A r }, for any i ^ the vectors A\x, . . . , A r x do not 
vanish simultaneously, and then p~ > as well as r] n > p~ > 0. 



i n 



3.5 Sets u> n and Q n 

Define for each n = 0, 1, 2, . . . the sets 



-Vl = |x € ^ m : P \\ x \\n = m & x ll^i x lln } i 

fi n = jx <E M m : /9 + ||x||* = max ||^ix||* | . 

By (|22h w n and fi n are the sets on which the value 

max^ ||^4.2*^ Hn 

\\ x \\n 

attains its minimum and maximum respectively. 
Lemma 3.8. The following relations hold: 

IMI'+l = {p n + (1 - PnjV^P") \\x\\n f° r X G W «, 



(24) 



X 



n+l 



/in + (1 - Ain)f?n /° + ) INI* / or ^ ^ ^n- 



Proof. The statement of the lemma is obvious for x = therefore in what follows 
it will be supposed that x G w n , x / 0. In this case (i24|l and the inequalities 
p~ < p + imply maxj ||^4jx||* = /9 _ ||x||*. From here by the definition (|23[) of the 
norm || • || n+1 we obtain 

ll^lln+l = Mnll^Hn + (1 ~~ Pn)Vn ^ max ll^-i^lln = (Pn + (1 — pn)Vn P ) ll^lln- 

For x G uj n the required equality is proved. For x £ Q n the required equality can 
be proved similarly. □ 

Lemma 3.9. For each n = 0, 1,2, . . . i/ie inclusions w n +i C w n , fi n +i f= /io/d. 

Proof. Let x 6 w n +i- If x = then clearly x £ oj n . Therefore in what follows it 
suffices to suppose that x ^ 0. In this case, by definition of the set u n+ \, 

max||^x||' +1 = p~||x||* +1 = p~ (p n \\x\\' n + (1 - Pn)^' 1 max||Aix|| n ) . (25) 



On the other hand by substituting || • ||* for the norm || • || n in (fl~8j) - (f2"0|) , and p , 
p n and r\ n for the parameters p~, A„ and 7„ respectively, we obtain the following 



estimate for max, ||^4jx||* +1 : 



max 

i 



\Aix\\l +1 > p n max + i 1 ~ Pn)v n V m ax ||Ajx||* . (26) 



Since by Lemma 13.71 p n > A > 0, from ([25]) . (f26|) it follows that p ||x||" > 
maxj ||^4jx||* or, what is the same, 

max^ 1 1 A-iX 1 1 ^ 

P — 1 1 „„ 1 1 • 



This last inequality by definition of the number p holds only for the elements 
x £ uj n . So, the inclusion cj n +i C uj n is proved. 



1 1 



Proof of the inclusion Q n +\ C Q n can be provided similarly, nevertheless for the 
sake of completeness prove it too. 

Let x G £l n +i- If x = then clearly x G Q n . So, consider further the case when 
x ^ 0. In this case by definition of the set fi n +i) 



max||Aix||' +1 = P + |M|* +1 = P [Hn\\x\\n + (1 - Pn)Vn m & x \\ A i x \\n • ( 27 ) 



On the other hand by substituting || • ||* for the norm || • ||„ in (fTHj) . (fT9|) . (|2TT) . 

and p~, fj, n and % for the parameters p~, X n and 7„ respectively, we obtain the 
following estimate for max, 

max ||^4jx||* +1 < p n max + (1 — p n )r]~ l p + max (28) 

i i " " i 

Since by Lemma I37F1 p n > A~ > 0, we see that (|27l) . (f28l) imply p + ||x||* < 
maxj \\A i x ||* or, what is the same, 

I maxj ||^4jx|L 
P ^ • 



By definition of the number p the last inequality holds only for the elements 
x G f2 n . Thus, the inclusion f2„+i C O n is also proved. □ 

Corollary 3.10. u = n n > w n ^ and = n n > r2„ ^ 0. 

Proof. By Lemma 13.91 {w n } is a family of embedded closed non-zero conic sets. 
Then the intersection oj of these sets is also a closed non-zero conic set. The same 
is valid for the sets {Q n }- D 



3.6 Completion of the proof of Assertion A2 

Choose non-zero vectors g G n n >oco n , h G n n >of2 n which exist by Corollary 13.101 
Then by Lemma 13.91 for each n > the following equalities hold: 

IMIn+1 = (/«n + (1 - Vn)rin l p~) IMI* , 

ii^ii;+i = (^ + (i-Mn)^ 1 p + ) \\h\\' n , 

From here 

IMI*=£nlMlo ; INI* = 6tWloi n ^ ' 

where 

n n 

c = n +0-- ^k)ri k i p-} , c = n k + ( x - A 1 *)^ v} • 

fc=0 fc=0 

The eccentricities of the norms || • ||* are uniformly bounded with respect to 
some Barabanov norm || • ||* (this fact can be proved by verbatim repetition of the 
analogous proof for the norms 1 1 • 1 1 n ) . Therefore the norms 1 1 • 1 1 ^ form a family, 
uniformly bounded and equicontinuous with respect to the Barabanov norm || • ||*: 

3 5 ± e{0,oo): S~\\x\\* < \\x\\l < 5 + \\x\\*, n = 0,1,2,.... 



1 n 



Then the sequences {||g||*} and {||/i||* } are uniformly bounded and uniformly sep- 
arated from zero, and the same holds for the sequences {£,„} and }. Let us 
show that the latter can be valid only under the condition p~ = p + . 
Note first that the inclusions rj k G [p~,p + ], valid by (|22p for all k, imply 



H k + (1 - n k )r] k l p < 1, k > 0, 
p k + (1 - Pk)v k 1 P + ^ i> ^ > o. 



(29) 
(30) 



If we additionally suppose that p < p + then the inclusions /U n € [A , A + ] and 
Vk G [p~ i P + ]i valid for all /c, will imply stronger estimates: 



p, k + (1 - n k )rj h 1 p < 

A+ + (l 



A + ) 



2p- 



< 1 if 7] k G 



(31) 



and 



A- + (l-A~) 



2^ 



> 1 if r\ k G 



• (32) 



Now, note that under the condition /> < p + infinitely many of numbers n k get 



into one of the intervals 



P 



or 



many indices k the estimates (|3ip are va 



/ "+p 4 



Therefore either for infinitely 

id while for the rest of them the estimates 
(|29p hold or for infinitely many indices k the estimates (|32p are valid while for the 
rest of them the estimates (|30p hold. Then in the first case £~ — >■ while in the 
second case — >• oo. 

Thus, in any case the assumption p~ < p + leads to the conclusion that the 
sequences and cannot be uniformly bounded and uniformly separated 

from zero simultaneously 

So, the proof of the equality p~ = p + is completed, and hence the iteration 
procedure (|5|), © is convergent. 



Max-relaxation iteration scheme 



In [2l|], for the same purposes, it was introduced the so-called max-relaxation proce- 
dure. We describe it shortly. Let 7(t, s), t, s > 0, be a continuous function satisfying 

7(t, t) = t, min{i, s} < ^(t, s) < max{t, s} for t ^ s. 



In [211 ] such a function is called an averaging function. Examples for averaging 
functions are: 

7(M) = ^A j(t,s) = Vts, j(t,s) 



t + s 



Given some averaging function 7(-, •), construct recursively the norms || • || n and 
• ||°, n = 1,2,..., in accordance with the following rules: 



1 o 



MRi: if the florin || • \\ n has been (ilvcddy defined compute the quantities 

I lUcLXj H^l^X 1 1 fi IXLciXj ||v4^x||t2 / i , / \ 

Pn = max |f— 7| > Pn = nun rr— ji , 7n = TIPn , Pn J! (33) 

MR2: define the norms \\ ■ \\ n +i and \\ ■ || n+1 : 

||x|| n+ i = max I \\x\\ n , 7 7 7 1 max ||Ajx[| n j , (34) 

\\ x \\n+l = lklln+l/||e|| n +l- (35) 

The max-relaxation procedure (f33|) (f35]) (the MR-procedure) possesses the same 
convergence properties as the LR- procedure [2ll ]. 



5. Examples and concluding remarks 

Several dozen numerical tests with 2x2 matrices were carried out with the help of 
MATLAB. Two of them, quite typical, are presented below. In the LR-procedure 
the relaxation parameter A n was chosen to be identically equal to 0.3, while the 
averaging function in the MR-procedure was taken as follows: A(i, s) = (t + s)/2. 

Example 5.1 Consider the family £/ = {yli,^} of 2 x 2 matrices 

^={1 0' A2= {-\ ?)■ 

The functions 3>i (</>), Hi (99), R n (ip), Rn((f) were chosen to be piecewise linear with 
3000 nodes uniformly distributed over the interval [—it, vr]. It was needed 21 steps of 
the LR-procedure and 22 steps of the MR-procedure to compute the joint spectral 
radius p{s^) with the absolute accuracy 10~ 3 . The computed value of the joint 
spectral radius is p(^) = 1.389. The computed unit sphere of the Barabanov 
norm || • ||* is plotted on Fig. [T]on the left. 

Example 5.2 Consider the family = {Ai,^} of 2 x 2 matrices 

_ ( 15/17 -16/17\ / 4/5 3/5\ 

Al "V4/17 15/17 j ' A2 ~ V-3/5 4/5 J ' 

Here the functions &i((p), Hi((p), R n (ip), i?* ((/?) were also chosen to be piecewise 
linear with 3000 nodes uniformly distributed over the interval [— vr, tt] . It was needed 
31 steps of the LR-procedure and 25 steps of the MR-procedure to compute the 
joint spectral radius p(s^) with the absolute accuracy 10~ 3 . The computed value 
of the joint spectral radius is p(s^) = 1.192. The computed unit sphere of the 
Barabanov norm || • ||* is plotted on Fig. [T]on the right. 

As is seen from these examples the computational "quality" of the above iteration 
procedures is approximately the same. At the same time similar steps in their proofs 
require different efforts and potentially may have different theoretical extensions, 
and now we are unable to predict which of these two algorithms might be more 
useful in the future. 

In conclusion note that the above algorithms allow to calculate the joint spectral 
radius of a finite matrix family with any required accuracy and to evaluate a pos- 
teriori the computational error. At the same time the question about the accuracy 



1 A 



Figure 1. Examples of computation of Barabanov norms for a pair of 2 X 2 matrices. 



of approximation of the Barabanov norm || • ||* by the norms || • || ra is open. It 
seems, the difficulty in answering this question is caused by the fact that in general 
the Barabanov norms for a matrix family are determined ambiguously. Namely to 
overcome this difficulty we preferred to consider relaxation algorithms instead of 
direct ones. Moreover, if to set A n = in (|6|) then, as demonstrate numerical tests, 
the obtained direct computational analog of the LR-procedure may turn out to be 
non-convergent . 

The question about the rate of convergence of the sequences {p^} and {Pn} to 
the joint spectral radius is also open. 
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